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^3 , The evolution of gauge invariant second-order scalar perturbations in a general single field infla- 

^^ ■ tionary scenario are presented. Different second order gauge invariant expressions for the curvature 

Ch ' are considered. We evaluate perturbatively one of these second order curvature fiuctuations and a 

second order gauge invariant scalar field fiuctuation during the slow-roll stage of a massive chaotic 
infiationary scenario, taking into account the deviation from a pure de Sitter evolution and con- 
sidering only the contribution of super-Hubble perturbations in mode-mode coupling. The spectra 
resulting from their contribution to the second order quantum correlation function are nearly scale- 
invariant, with additional logarithmic corrections with respect to the first order spectrum. For all 
scales of interest the amplitude of these spectra depends on the total number of e-folds. We find, on 
-y~, I comparing first and second order perturbation results, an upper limit to the total number of e-folds 

^ . beyond which the two orders are comparable. 
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(3 ■ I. INTRODUCTION 

^^ ' One of the observational arguments supporting inflation is the evidence of nearly scale- invariant primordial adiabatic 

O . fluctuations with a mostly gaussian spectrum, as shown by COBE [l| and WMAP [1,|3- In its simplest realizations, 
O^' adiabatic cosmological fluctuations are simply related to linear fluctuations of the inflaton generated during the 
Vh . accelerated stage: such fluctuations have a gaussian spectrum. 

Ol)' Non-linear fluctuations inevitably carry a x-squared distribution, which may contribute to non-gaussianities in the 

^ I 3-D power spectrum and in CMB anisotropics depending on the amplitude and spectrum of these non-linearities (see 

: '~j ■ [Jl for a review) . A study of fluctuations beyond the linear level is therefore interesting not only for the issue of the 

KS ' stability of inflationary dynamics, but also because of its connection with the predictions of non-gaussianities in the 

2 . matter power spectrum and in CMB anisotropics. Even in absence of non-linearity in the potential for the scalar field, 

■ - - ' the generality of intrinsic non-gaussianities is guaranteed by the non-linearity of the Einstein equations, a fact which 

has motivated a great interest both from the theoretical and observational point of view [j| . Our choice of calculating 

the intrinsic non-gaussianities for a free massive inflationary model is therefore a case of primary interest, since in 

such a case the whole non-linear sector is given by gravity. 

In this paper we shall focus on the perturbative evualuation of second order gauge-invariant (GI henceforth) 
curvature fluctuations. For this purpose wc discuss different second order GI measures of curvature perturbations 
(see also |a,lyjl3jBBIl3ll3)i which to first order coincide with the curvature on fiat slicing hypersurfaces. We then 
compute the correlation function for differents second order GI variables, finding, for all scales of interest, a growth 
in the spectrum amplitude as the total number of e-folds which infiation lasted increases. A dependence on the total 
number of e-folds was also found for the production of test scalar fields with a mass smaller than that of the infiaton 
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The paper is organized as follows. We present the formalism and different gauge-invariant measures of curvature 
perturbations in section II. In sections III and IV we evaluate the quantum correlators for different second order GI 
variables with both analytical and numerical methods. We conclude in section V and in appendices A e B we add 
useful technical formulae for gauge transformations and for the analytical results given in section III. 

II. GENERAL CURVATURE PERTURBATION 

Let us consider inflation in a flat universe driven by a classical minimally coupled scalar field with a potential V{4>) 
and follow the notation in Il3l|. The action is given by: 



S = / d'^xC 



Yla-yd,^9.^^vi^) 



(1) 



where £ is the Lagrangian density and we shall study the chaotic model V{(J3) = ^4''^ [3 (see [T3 for a second 
order calculation in a different inflationary model). We shall consider, up to second order in the fluctuations, both 
the inflaton: 

(/)(t, x) = (j){t) + ip{t, x) + (^(2) (i^ x) 

and the metric of a flat universe (ds^ = gf^udx^dx'^): 

goo = -1 -2a-2a(^^ 
90. = -i(A.+/3?) 



,2 



5,j - 26,, (^ + 7/.(2) ) + D,, (^ + ^(2) ^ ] ^ (2) 



where Dij = didj — 1/3 V^(5.y . As is clear from the above we limit ourselves to the scalar sector of metric perturbations, 
neglecting vector and tensor perturbations. 

In the above formulation the gauge is not fixed, and one can eliminate two scalar functions among the four metric 
coefficients and the inflaton fluctuation. We wish to extend to second order the GI potential for the curvature 
perturbation. For this purpose, let us note that to first order the so called curvature (potential) perturbation is given 
by 

^ = ^ + lv2£; (3) 

and th is q uantity is not GI (for the general gauge trasformation see appendix A). We now consider the comoving 
slicing [l^j which is defined to be the slicing orthogonal to the worldlines of comoving observers: they are free- falling 
and the expansion defined by them is isotropic. That is the observers do not measure any flux of energy {T ^ =0) 
and for the universe under consideration this corresponds to having ip = 0. The transformation from a general slicing 
to a comoving slicing with tp = is: 

V'coni = ^- e"i)0 = ^ e"i) = -r (4) 

and one obtains 

R^''>^i^,orn^i'+J^+l^''E, (5) 

which is the first order GI comoving curvature (potential) perturbation. 

To second order the situation is more involved due to the presence of terms quadratic in the first order. We are 
interested in the definition of GI quantities including second order. Their construction is not unique and we shall 
therefore construct three different GI second order curvature (potential) perturbations. 



A. The comoving curvature potential perturbation Ra 

This variable is obtained by repeating consistently up to second order the procedure used for the first order comoving 
curvature (potential) perturbation. Starting from 



^(2) = ^(2) + iv2s{2) , 

6 



(6) 



to second order the comoving slicing is caracterized hy T ^ =0 and T ^ = 0, that is <^ = and (^^^^ =0. In 
this case, in order to fix all the degrees of freedom, we restrict ourselves to considering only an infinitesimal time 
trasformation up to second order, which implies taking e/-^-, and €%-, equal to zero. From the first order we obtain the 
condition in Eq.l^J, while from the second order we have 



'pioL = V 



(2) 



£(!)¥'■ 
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^?i)(^(i)^)' 



t(2)^ 
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^'^'-4)^ + ra)(^a)'^)' 



(7) 



Thus one obtains 
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6a d, 



d'Pd.ip - -V^ (Pip) 



Id'd^ 



4 V2 



-2H%D^jE 



(h a^ d)'' la (h 



(8) 



Let us note that this quantity is analogous to the second order curvature perturbation on a uniform density hyper- 
surface defined in Q- These two quantities agree (apart from a sign) in the long wavelength fimit (see also [all3)- 



B. The curvature potential perturbation Rb 

We give the second order GI curvature perturbation potential obtained from a form analogous to Eq. (jSJ (the first 
three terms of the following equation), and add a minimal set of terms necessary to obtain a second order GI variable. 
Using the general transformations of Appendix A, we obtain 

42) ^ ^(2) ^ H^i2) ^ 1^2^(2) 

S 6 



— 2H^ + H - H^ 

'6^ \r +2(V^ 



t/j + 2Htlj - 



H 



-ip 



24 



d^(3d'P 



-i/3Aj/3 + 2i'D,,E + ^^D,jE 



(9) 



This quantity is the true GI generalization of T?.'-^-' computed in [g|, which also takes into account infinitesimal spatial 
(scalar) transformation in addition to the infinitesimal time translations considered in [3. 

C. The comoving intrinsic curvature potential perturbation Re 

In order to introduce the third curvature (potential) perturbation, let us first consider the intrinsic 3-curvature ^^' R 
associated with a foliation of space-time with constant proper time t for a fiat Universe. Using cq.(49) and cq.(58) of 
la we obtain: 



(3)i?, 



1 



2d'd,Cl ~ 2d'd,C] + 4C'^ {-2djdkC'l + dkd'^C + d,d,Cl) ~ {20^0^ ~ dfil) {2d,C'' - d'Cf) 



+d*'C'^{?,duC,j-2djC,k)\ 



(10) 



with 



a 



ij 



-5,, (V' + ^(')) + iA,(i? + i?^'^ 



In our notations the intrinsic curvature is: 

(3)^ = (3)i?(l) +(3) i?(2) 



+ ^5» (i' + l^'E^ d. [i; + \v'E^ \v'EV^ (^0 + f^'E 
-^d'd^Ed.d, U' + I^^e) - ia' U + 1^^e\ d,\/^E - ^d'\/^Ed^\/^E 
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d'd^d^Ed^djdkE 



(11) 



Instead of using Eq. lO, we now consider the potential ijjfj of the intrinsic curvature to second order ^'^'R^'^\ given 

by (^)ii;(^) — ^V^f/'^ . We then obtain the comoving GI expression by going from a general slicing to a comoving 
slicing: 
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(12) 



D. Comparison among curvature perturbations in the long- wavelength limit 

The preceding three expressions undergo a considerable simplification in the long- wavelength limit. For large scales 
we obtain: 
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,<;>-_,..,« .>^U,..H,)-^.,,t,-n'-"-,!it 



H H > 



(13) 



(h 2 \ d> I \ (b 



(14) 



R. 



(2) 
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-H^- — + — 4- 



^i*-^f rn*-f-Mn*^f->- *+f 



v' 
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in this limit the difference between R''^' and i?^'' is given by [3: 
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2H^ + H -H^ 
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Let us note that, as for the first order R^^\ R)^ also is constant in time in this limit 

Further on using different parametrizations of the metric many other gauge invariant quantities can be constructed. 
For example from |ll| one can define another variable -R^^, which is related to i?^ on large scales by 



R 



(2) 



R 



(2) 
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i?(i; 



(17) 



(2) 

We also consider R^ since it is directly related to a geometrical quantity: the intrinsic 3-curvature. 



III. LONG- WAVELENGTH SECOND ORDER CURVATURE PERTURBATIONS IN UCG 



As in one of our previous papers on back-reaction [l^j, we choose to work in the uniform curvature gauge (UCG), 
which can be generalized straightforwardly to second order through the conditions: tp = tp^^^ = and E = E^'^'^ = 0. 
In this gauge, the evolution equation for inflaton fluctuations is regular also during the coherent oscillation of the scalar 
field [13. In the UCG and in the long wavclcnght limit we have, for the given second order curvature perturbation 
potentials, the following results: 



^'^— .-'-'- T.^--^-' 
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ij2 



-ip 



(19) 



i?: 



(2) 



H 



-LP 



(2) 



H 



-wip- 



H' 







3H' 



d'lpdiip 



(20) 



On subtracting the average of the quadratic piece (in order to also have a zero average for the non- linear piece) , it 
is useful to introduce the parameter of nonlinearity Jnl jllj : 



R^Rl + ^/nl [Rl - {RD] 



(21) 



where i?L is the Gaussian part of R (that is the first order part R^^^). Using the first order comoving curvature 
perturbation together with the first two second order definitions given in Eqs. H18|l . H19|l . we obtain for a ^0^ 
chaotic inflation and using the approximation 99 = —-^i/? and Ea. H47|) . to leading order in the slow- roll parameter 
e = —jU, the following results for f[fi^ : 



rA 

/nl 



/nl 
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Let us note that we cannot evaluate /^^ ii^ analogy with |lll | because the non-local terms are not subleading in the 
long-wavelength approximation. 

In the next two sections we proceed with the evaluation of the spectrum of the non-linear corrections to the curvature 
Ra and of the second order GI field fiuctuation Q'^^ defined in [l^- The evaluation of the quantum correlator of 



second-order GI variables (see also [l9|) involves a sum over momenta, which is plagued by ultraviolet divergencies. 
The adiabatic subtraction used in the calculation of the finite part of energy and pressure of fluctuations [l^, U^ is 
roughly equivalent to only considering super-Hubble fluctuations. Motivated by this analogy we shall also use this 
criterium in the present case: all the integrals which follow have the physical Hubble radius as an upper limit of 
integration. The same integrals also need a lower limit of integration which defines the region of nearly scale invariant 
(red tilted) spectra generated by infiation: we choose this lower limit of integration, a comoving momentum i, of 
order of the Hubble rate when infiation begins, Hq, as previously found in analytical [20i | and numerical investigation 
[l^ (see also |21|). At the present state of our knowledge (i.e. without any theoretical input for the choice of £), we 
believe that this is a natural choice as compared to considering £ of the order of the present comoving Hubble radius 
as is done in 1^. From now on wc shall restrict ourselves to the slow- roll stage of ^(ff^ chaotic inflation, during 
which H ~ —rv? ji holds. 



A. Analytic Evaluation of curvature IVj^ 

Within the above approximation from cq. H18|) , to leading order in the slow-roll parameter and using Eq. (28) of our 
[la. one obtains 



^V'^') 



si?«^ 



(22) 



and 



rI^ ^R^^'^^\-\ + 0{e)\ . 



(23) 



From the above equation and Eq. (|17|l it is evident that Rlr is of order a slow-roll parameter times the square of the 
first-order perturbation. 

Let us consider the quantum correlation function for this second order gauge invariant variable 



(0|i?(f)(x)4^)(y)|0) = (l + O(6))(0|(i?W(x))'(i?W(y))'|0) 



:^(0|(^(i,x))^(^(i,y))^| 



(24) 



The quantized scalar field variable ip is given by 

(^(t,x) = -n TTTTTT / rfk 



Mt)e^- a^ + pl{t)e-^'^- a\ 



(27r)3/2 
in the slow-roll approximation and on large scales (but < 2ti / £) we consider [ic 

•^k{t) - 



(25) 



(26) 



with 



iJ(tfc)~i7o l-2eoln 



Hn 



1/2 



(27) 



which is the Hubble parameter when the fluctuations crosses the horizon and eo = e(t = 0). The power spectrum of 
first order curvature perturbation, at large scales and in the UCG, is therefore: 



P^a,(A:)^^|i?, 



(1)|2 






(28) 



It is important to note that by using the expression 127|) for H{t]^) we obtain a power spectrum of curvature per- 
turbations in H28|l having the correct value to leading order of the spectral index and exhibiting running during 
slow-roll. 



(2) 

The physical content of the correlation function for R)^ is conveniently described by the power spectrum associated 



i?i')(x)4^)(y)|0) = j,;^ / d^ke^'^-^^-y^lRf^it)]^ (29) 



with its fourier transform 

1 

which becomes, using Eq.(j25 

eiaH - fc2)e(fci ~ O0(fc2 - exp^('^^+''^)(''-y) + / d^k^d^k2^^^^^^^^^ 

e{aH - ki)Q{aH - k2)Q{ki - l)Q{k2 - I)} , (30) 

where we have considered only modes which are super-Hubble, but greater then comoving infrared cut-off Z, using the 
Heaviside functions O. We can perform a Fourier transform with respect to r = x — y of Ea. (|30|l and obtain 

e(fc - Oe(|P - k| - /) + (2^)3^(3) (p)Qo (31) 

where 

Qo = ^^^ I d^fcid3fc2:^^£(|a)!_l_e(„ij _ k,)e{aH - k2)Q{k, - l)e{k-2 - I) . (32) 

For p ^ and on evaluating the integral in Ea. H31|) to leading order and in the momentum window given by the 
condition I << p << a{t)H{t) one obtains the expression: 

which is time independent in our approximation. The coefficients Ai are given in appendix B, g (-] = C ( (in ?) - 

""(..fe)-'5((>«^^^)*(.(7fel)> 

Substituting for the coefficients Ai we obtain the following result: 

l4'il' - ^TTT^ (-2.8405iiJ(ip)6 + 19.3975H{tp)* - Um41HH{tp)' + 9.72737J2 + 4^H{tp)' 

,., I , ,„,y. (,., I) % «,„). (,„ |)% ,, (1) , , (_E_) } . ,3., 

On taking only the leading term for the coefficents of the different powers in ln(p/Z), we obtain: 



^P 327r2 Ai^j F I ™2 m* I \ U 3 \ U 

1 , ,,^„ f H(uy H(t„Y f p\ m^ ( p\^ 16 to2 



where in the last equality we have used Eq.jSSJ. This last expression shows compactly one of the main results of our 
investigation: the spectrum of this second order GI curvature perturbation is proportional to the first order spectrum 
through logarithmic corrections which encode the scale I at which inflation started. Without these logarithmic 
corrections the second order curvature perturbation would be completely local and much smaller than first order 
terms since for the scales of interest H{tk) « Mpi. 



B. Analytic Evaluation of Q^^^ 

Let us consider the second order GI scalar field fluctuation Q^^^ (which is a possible second order generalization 
of the first order Mukhanov variable [23) defined in [l^ as the second order field fluctuations on uniform curvature 
hypersurfaces in the longwavelength limit. The uniform curvature hypersurfaces is defined as the slicing with ip = 
i/j^^^ = 0, and one obtains the following result |7l[ia|: 



)(2) = ^(2) + ^^(2) 



H 






(36) 



In the UCG; as for the first order, this GI variable is simply the second order field fiuctuations (p'-^\ studied in our 
previous paper ^^. Thus the equation of motion for this second order GI field fluctuation in the UCG, to leading 
order in the slow- roll parameter and in the long- wavelength approximation is given by [l^j : 



i(2) 



i(2) ^ S77^n(2) = ,7^2 



Q^^' +3HQ^^> +3HQ 






(37) 



It is interesting to comment on the self-interaction mediated by gravity corrected for the inflaton. Just as the 
feedback of metric perturbation vanishes in first order perturbation theory for ^ = (see Eq. (7) of [13)7 it also 
vanishes for the self-interaction to second order in Eq. \i7\ . 

Integrating the above in two steps we obtain: 



)(2) 



Q^^> + 3HQ^^' = m^ / dt 



,, yft^x)^ 
2H{t') Ml 



(38) 



and 



)(2) 



(i,x) 



a{tf 



'„W\3 



dt'a{t 



dt" 



" -^^2 



2H{t") 



,(^(t^_x 

M?. 



(39) 



Let us consider the quantum correlation function of this second order Q*-^' g^-uge invariant variable 
(0|(5'^^(i,x)(5(^'(i,y)|0). The quantized scalar variable Q^'^^ is given by 



)(2) 



(t,x) 



(27r)3/2 



dk 



Qf{t)e^--b,. + Q)^^*{t)e-^'^--bl 



(40) 



The physical content of the correlation function for Q^^-* is better described by the power spectrum associated with 
its fourier transform 



(0|Q(2)(t,x)Q(2)(t,y)|0) 



1 



(27r)3 



d3fce^k.(x-y)|Q(2)(^)|2 



and using Ea. l39|l and Ea. l2()|l . one obtains 



(41) 



)(2) 



1 I 4> ra 



(0|Q(^Ht,x)Q(^)(t,y)|0) = — -^-^ 



a{tf \2Ml {2nr 



d^kid^k2 



ki ki 2 



'„f^'\3 



dt'a{t') 



1 



^^"i?(FF®^''^*"^^^^"^ fci)e(a(t")iJ(i") - fc2) 

d:'kid^k2 



e(fci - l)Q{k2 - /)exp*(''i+''2)(^-y) 



-3,. ^3;, H{tk,f H{tk,f I 

kf fc3 4 



dt'a{t'y / dt 



1 



dt'aitr I dt"j^^Qiait")H{t")-k2) 



H{t"f 
Q{ki-l)Q{k2-l) 



Q{a{t")H{t") - fci) 



(42) 



On performing the Fourier transformation w.r.t. r = x — y of Ea. H42() one finds 



(j) m? 



IQ^'WI = 777W iiTT 



a{tf \2Ml (2^)3 



fc3 |p - k|3 2 



df'^ii'f I dt"——^e{ait")H{t")-k) 



H{t'' 



Q{a{t")H{t") - |p - k|)]' e(fc - Oe(|p - k| - + (2^)3<5(3)(p)Qo 



(43) 



where 



1 



^"-;^Ui| (2^)^^^^^^^^^ 



H{tk,)^H{tk,fl 



kl 4 



dt'a(t')''' / dt"—^Q{a{t")H{t") - fci) 



i?(t 



1 



ds'ais'f / ds"-^^^-^e(a(.")i?(s") - fc2) 



e(fci - l)Q{k2 - I) 



(44) 



On considering p ^ Q and replacing the O function by the domain of integration, it is possible to perform the time 
integrals exactly: 



I I Mt)H{t) 
2 Ml) {2^2Ji 



\Q^:\t)? = l^^l IT^Al dkl dyQ{a{t)H{t)^{p' + e-2kpyf'^-)Q{{p' + e-2kpyf'^^-l) 
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H 



k (p2 + A:2 - 2fcpy)3/2 

1 _3iim!. f 1 

-e 2 H 



"^'^W'TfA '+^w'° 



i? , (p2 + fc2 - 2kpy) 



1/2 



367J2 
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2H 



\f{k,p,yY \ 2H 



1/2 



^,1 3^ 
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2' 2 ij 



2' 2 i/ 
1 if{k,P,y) 



Ho J 

1 _3/(fc^p^2^^ 
2' 2 ij 



2\ 1 



2' 2 H 



(45) 



where f{k,p,y) is given by 

f{k,P,y) = Mm 



"" '+4'°i 



1/2 



, i^o I 1 + 2-2 In 



F {p^ + P - 2fcp2/)i/2 



1/2- 



It is also possible to obtain an explicit approximate espression for the integral starting from Ea. H43() subject to the 
condition / << p « a{t)H{t) and the approximation I « p « a{t")H{t"). In particular we find 



1 „. 1 



>^^^W.^(|^jii.o^^|EA,(lnf 
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(46) 



where, as in the previous subsection, the coefficients Ai are discussed in the appendix, .9 ( - ) = ^ [ (In f ) ~ ) ^-nd 
Let us also investigate a much more crude approximation, obtained on neglecting Q^'^' in Eg. 1)38(1 and considering 



the infrared limit 93 = —^f, which leads to the expression 



Q(2) ^ L 



1 6 1 



4 H M2 
pi 
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^£^(1)2 
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(47) 
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The correlation function for Q?' turn simplifies to 

2^e(aH - h)e{aH - k^)Q{k, l)Q{k2 I) exp'(''^+''^)(--'') + j d^k.d'k^^^^^^^^, 

e{aH - ki)Q{aH - k2)e{ki - l)Q{k2 - 0} • (48) 



So, on proceeding as the previous subsection, for p 7^ 0, to the leading order and in an "approximate" fashion with 
the condition / << p « a{t)H{t) one obtains 



I \ I H^ I \^ . L W fl 



ie?w= 5l^ s^^? i:A(.nf)V4y.H^(oW) ■ <-' 



Substituting the coefficients Ai we obtain the following result: 
2 



Upit)? ^[j^^] ^^'^ |-2.8405-liJ(ip)« + 19.3975H(tp)4 - umAlHH{tpf + 9.7273i72 + 4-^iJ(i,)8 



On taking only the leading term for the coefficcnts of the different powers in ln(p/Z), we obtain: 



lo?'(Oi' . _L.,,,,,P ...,.S|Z. 3.fM (,„.), „^ (,„.)% -i|(,^ ,., 



Analogous considerations to those given after Ea. (|35|l are also valid in this case. Also, as we can see from Eq. (|50|l . 
one obtains a spectrum which, according to Eqs. ^2'6\ and (|47|l . slowly increases in time during slow-roll as e" , whereas 

Rji/ is constant. Again, we note that both the GI second order field fluctuation Q^^^ and the second order curvature 

R\ have a nearly scale invariant spectrum with a logarithmic corrections. 

IV. NUMERICAL ANALYSIS 

Let us introduce a few parameters which are useful for the numerical analysis. The first is Ntot = In (o(i/)/a(ii)), 
the total number of e-folds, thus we need to specify the time associated with the end of inflation. This definition may 
be obtained on requiring the Hubble parameter H{t) to be equal, for example, to the inflaton mass to. This happens, 
in our chaotic model for the time tf when the slow-roll parameter e is equal to 1/3. We then solve the equation 



obtaining 



e(i) = e (52) 



(53) 



i=- 


-f-; 


1 
m 




Ho 


1 3 


h = 


H \ 


l-H 


by 






Ntot 


iHl 
2to2 


3 

2' 



and the final time, defined by e = 1/3, is 

(54) 
n \ n 

corresponding to a total number of e-folds given by 

(55) 
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It is also useful to define the variable N{t), which is equal to the number of e- folds from the end of inflation, 



a{tf) _ 2,H{tY 3 



N{t)=\n^ = -^^--. (56) 

In our analysis we are interested in the fluctuations that cross the horizon at A^* — 55 e-folds from the end of inflation. 
The time i* associated with N^ = In [a[tf)/a{t^)) = 55 is: 



H \ H H 
and the comoving wave-number which crosses the Hubble radius at A'^* is 



u = ^^- J -I In. (57) 



k, = a(i,)iJ(i,) = e~^-^-^' \J-2,H - 2HN. . (58) 

If we wish to take a "photograph" of a fluctuation with k. on super-Hubble scales during inflation we should take: 

->e> i . (59) 

3 27V* + 3 ^ ^ 

We now proceed with the numerical analysis taking m = 10^^ Mpi. In particular we are interested in the comparison 
of second order fluctuations with first order ones, in order to understand for which condition first order perturbation 
theory breaks down. We shall also plot the shape, in k, of the spectrum of fluctuations. 

A. Second order GI scalar field fluctuation Q^'^' 

In figure 1 we consider, for two different values of the initial Hubble parameter Hq {7m and 14771, respectively), the 
logarithm of the second order contribution to the power spectrum of Q"^: 

as given by Ea. (|50|l . Ea. H46() and the numerical solution of Ea. (|45() . in order to see the improvement of the roughest 
analytical approximation, with respect to the exact numerical calculation with increasing Hq, and the shape of the 
spectrum. We give the figure at an instant during inflation for which e = 1/10. 
Corresponding to those two values of Hq we have the following data 

• For Hq = 7m we have iVtot = 72, tf = 18/m, t. = 2.58805/m and k./m ^ 1.48247 10^. The time at which we 
"photograph" the inflation is t = 15.5228/m. 

• For Hq = Um we have Ntot = 292.5, tf = 39/m, U = 23.588/to and fc,/m = 8.56876 10^"^. The time at which 
we "photograph" the inflation is t = 36.5228/m. 

In order to compare the F^'^\k,t) given by the numerical solution of Eq.|03J) with the first order contribution to 
the power spectrum 

given by Eq.J^H), we give two different types of figures. In figure 2 we compare the values calculated at fc, as a 
function of the total number of e-folds Ntot, at the end of inflation t = tf. From the figure on the right (which differs 
from the figure on the left only because of the bigger range of Ntot) we see that second order effects are of the same 
order as first order effects for Ntot near 30000 which corresponds to a Hq near 141m. Therefore for such values first 
order perturbation theory breaks down 

In figure 3 we show the dependence of the fluctuations on the modes, crossing the horizon N e-folds before the end 
of inflation, as seen at the end of inflation. We consider three different cases corresponding to values of Hq/iti equal 
to 7, 14 and 30, compare F^^'>{k]\j,tf) and F^^'>{k]\j,tf), and vary TV from 55 to 1 with: 



tN 



'^-^-^\/-3H -2HN (60) 



which is the mode of the fluctuation that crosses the horizon at N e-folds before the end of inflation. As before the 
first and second order fluctuations are of the same order of magnitude for Hq = 141m. 
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login (fc/m) 




logy)(fc/m) 



FIG. 1: We plot for Ho = 7m (on the left) and Ho = 14m (on the right) the evolution with respect to logm (^) of F^'^\k,i) 
using the analytic approximations Eg. 150^ Mashed line) and Ea. l|46|l (dot-dashed line) and the numerical solution of Ea. H45|l 
(solid line). 



FW(fc.,t;) 



FW(K,tf) 



200 400 600 800 1000 1200 



Ntot 



-10 



-12 



-16 . 



F''\k„tf) 




Ntot 



FIG. 2: We plot F^^\k^,tf) and F^'^\k*,tf) with respect the total number of e-folds Ntot 



B. Comoving curvature perturbation R 



(2) 



Let us proceed in a analogous way for the second order curvature perturbation (potential) R\ . As before in figure 
4 we shall consider two different cases for inflation with Hq equal to 7ni and 14m and we shall plot 



^^'HA)=logro(|^l41r) 



the logarithm of a second order contribution to the power spectrum, as given by the Eq. (|34|) , to exhibit the shape of 
the spectrum. We give the picture an instant during inflation for which e = 1/10. 



F^^\kN,tf) 



-N 



-50 



-40 -30 



F'^^\kN,tf) — — — — — ^__ 




12 

14 

16 
\ 

1 



FIG. 3: We plot _F'^^(fcjv, i/) and F^^\kN,tf) for Ho = 7m (solid line), H, = 14m (dashed line) and H, = 30m (dot-dashed 
line) with respect A*'. 
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-14.6 



S^^\k) 




-12 


/^ 






12.5 




-13 




13.5 





logio(fc/"-i) 



20 40 60 80 100 120 log^^(k/m) 

FIG. 4: We plot, for Ho = 7m (on the left) and Hq = 14m (on the right), the evolution with respect log^Q (i) of 5'^^(fc). 



Ntot 



200 400 600 800 1000 1200 





FIG. 5: We plot 5''^'(fc.) and 5*^'(fc.) with respect to the total number of e- folds Ntot 



For the sake of comparing S^'^' (fc) with the first order contribution to the power spectrum 

5«(fc) = logio(|^|i?«|^ 

given by Ea. H28|l . we give two different kinds of figures. In figure 5 we compare those values calculated at A;*, as a 
function of the total number of e- folds Nt^t, at the end of inflation tf. From the figure on the right (which again differs 
from the figure on the left only because of the larger range of Ntot) we see that for this case the second order effect 
is of the same order as the first order effect for Ntot near 7560 which corresponds to a Hq near 71m. Therefore first 
order perturbation theory breaks down for a value of Hq smaller than that obtained before and we have a stronger 
limit. 

In figure 6 we again consider three different value of Ho/m, namely 7, 14 and 30, compare S'(^^(fcjv) and S''^-'(fcjv), 
and vary N from 55 to 1 at the end of infiation. As before the first and second orders will be of the same order of 
magnitude if we take Hq = 7lm. 



S^^\kM) 






-N 


— =5re- ^4iL___:^30__^ -20 


-10 










-"^o 




~ ~ - ^ _ 






S^^\kN) 


— — 


^-^ 


^2 
^-14 




"~^^^ 



FIG. 6: We plot S'^'(fc]v) and 5^^'(itjv) for Ha = 7m (solid line), Ho = 14m (dashed line) and Ho = 30m (dot-dashed line) 
with respect to A*'. 
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Finally wc wish to compare the critical initial values of Hq which we obtain in this paper, with two different 
limiting values. The first one is the value H^ obtained by requiring that back-reaction in the UCG gauge will 
become important by the end of inflation [i3jl2J|j when H ^ m. Such a value is: 



Hbr 



167r 



2^1/6 2/3 ,,.1/3 



m^'-'M. 



pi 



(61) 



The second one is the value obtained by the self-reproduction argument [2£ 



H., 



(Stt 



,3/4 



1/2 



1/2 71.1/2 



(62) 



This comparison is given in Table I. 



Ho limit for Rf 


71 m 


Ho limit for Q^^^ 


141 m 


Hhr 


108 m 


Hsr 


1449 m 



TABLE I: Comparison of the values of Ho for which first order perturbation theory break down with //br and Hs, 



We stress that neither of the two values of Hq or of iJbr are related to the self-reproduction scale [23 • The values 
of Hq obtained by considering the breakdown of linear perturbation theory in this paper are of the same order of 
magnitude as iJbr (sec also |2g for an analogous value). This fact is not surprising: both these calculations include 
second-order cosmological perturbations. The results for back-reaction were considered ambiguous since it is not 
simple to demonstrate the slow-down of the inflationary expansion in terms of a GI quantity. On the other hand in 
the present paper, the breakdown of first order perturbation theory is found for GI variables in the long-wavelength 
limit. 

The results summarized in the Table I and derived from figures 2, 3, 5 and 6 can also be obtained directly from 
Ea. (|35|l and Ea. (|51|l . In those cases and for Hq >> m on using Ea. (|60|l the leading logarithm values (that is the term 



with(lnf)')of|i?(f^j2 



and \Qp^J,itf)\'^ are given by: 



m 



(2) I 



\Q^^I(tf)\' 



/|^2 l^-^N 



H^ 



47r2 
3 



M2m4 



167r2 



<tff\vpAtf)\' 



H^ 



Mp2m4 



(63) 
(64) 



and the first and second orders are comparable for 



Ho 



2n Mpi 



1/3 



(65) 



and 



Hq 



/ Utt Mp 
\V3 m 



1/3 



(66) 



respectively, where Ea. (|65|l nearly give the same value as Tabic I (namely near 71m), while Ea. (|66|l only gives the 
same order of magnitude (actually 130m), which is not surprising considering the approximations made. 



CONCLUSIONS 



We have studied second order cosmological perturbations for a chaotic -^0^ inflationary model, and considered 
three second order GI measures of curvature perturbations. 

For the GI curvature perturbation studied and for the second order GI scalar field flcutuation presented in [l3| , we 
have found, for all scales of interest (see, for example, Ea. (|63|l and Eq.ljMJ), that the amplitude of the spectrum of 
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quadratic corrections grows with the total number of e- folds A^tot- For Q^^^ we also obtain a mild time dependence 
associated with the growth of the slow-roll parameter e. Since both the quantities studied here are GI, this dependence 
on A^tot cannot be a gauge artifact. 

The spectrum of those second order GI variables obtained by Fourier trasforming the curvature-curvature quantum 
correlation function is nearly scale- invariant, with additional logarithmic corrections with respect to the first order 
spectrum. On comparing the first and the second order contributions for the two different cases considered, one finds 
two different values for the initial Hubble parameter Hq beyond which first order perturbation theory breaks down. 
We found limits which are of the same order of magnitude as that one found by back-reaction in the UCG |13j . 
Neither the back-reaction limit nor those found here have anything to do with the self-reproduction scale. 

In the curvature-curvature quantum correlation function, the cross terms involving first and third order perturbation 
may also contribute. The third order curvature perturbation is clearly beyond the scope of the present paper and will 
be the subject of future work. Unless such a term cancels or dominates the growth we have found, our results imply 
that in single field models the amplitude of the spectrum of intrinsic non-gaussianities generated during the slow-roll 
stage will depend on the total number of e-folds which inflation has lasted. 
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VI. APPENDIX A: GAUGE TRANSFORMATIONS 



The gauge in our formulation (see Eq.(|5|) is not fixed, since one can eliminate two scalars among the metric and 
the inflaton degrees of freedom. After a general analysis we shall be interested in making a gauge choice in order to 
match previously studied situations. An infinitesimal coordinate trasformation to second order is |27l |: 



+'^ri)+ 2 Ki).'^*^('i) +*^r2)j ' 



(67) 



where e(i) and 6(2) are the coordinate changes to first an second order, respectively, and eL) 
It induces the following change in a geometric object T = T^^^ + T^^^ -I- T^^^ : 



'£(2) and e^^) 



-(I)- 



-(1) 



fW ^ j^(i) _ £^ 7^(0) 



(68) 



.(2) 



;^(2) ^ j.(2) _ ^^^^^y(i) ^ _ ^jr2^^^j.^ _ c^^^^To 
which leads to the following general gauge trasformation for our scalar perturbations. To first order we have: 



(69) 



v = </'-£(!)'; 



(70) 



a = a — e 



(1) 



/3 = /3--e?i)+2aefi 



V; = ^ + iJe^i) + -V^e^, 



(71) 
(72) 

(73) 



E = E -2e1 



(74) 



To second order one finds: 



^(2) 



V 



(2) 



ii)V 



,0 



_0 



(1) Vni)'' 



-(2) 



dwd'el) , 



(75) 



d(2) 



v(2) _ .0 



^-~^-^-\^ + \ivA^ + (!^ 
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(76) 



^ -^ --^(2)+«e(2)+2^^(e(l)) +^ 



-a,e^i)eOi) - eOi)9,/3 - He\^)d,l3 - e^.-^d,/] - -ad.ef,) 



(77) 



^ 



(2) 



?A(2) _ gO^j ( 2i7V' + V- 



Y«(2) Y«(l)'^(l) 



,0 2 



iJ + 2H^ 






vh-' 



(2) 



12 



(78) 



^^^^ =£^^^^ + :^/3^?i)-:i(^?i))^-<2) + >4i 



1 



(1) 



3 s^a^' 



2(V2) 






iHel^D^.ely 



(79) 



VII. APPENDIX B: Ai COEFFICIENT 



The coefficients Ai are given by: 



Ao(fc) = mel (C(3) - l)Ci(fc) - 10.198 eg C2(fc) + 24.656 62^3 (A:) + eq f g - |^2 j Ci(fc)C4(fc) 

/„2 \ 

-9.30635 eoC2(fc)C4(fc) +33.333 eoC'3(fc)C4(fc) - 



4 C2(fc)C4(fc)' +6.07344 C3(fc)C4(fc)' 



Ai(fc) = 4Ci(fc)C4(fc)2-8C2(fc)C4(fc)2 

A2{k) = 4eoCi(fc)C4(fc)-8eoC2(fc)C4(fc) + 2C2(fc)C4(fc)2+8C3(fc)C4(fc)2 

^3(fc) - 3 6gCi(fc)-ye2^2(fc) + 3C3(fc)C4(fc)2 



(80) 



wlicre 



Ci(fc) 






Hn^ 



i/n^ 



^0 -'-'0 

C2{k) = -2eo f^^ll^ - 260 



7^2 

-"0 



C3(fc) = 
C4(fc) = 



Hi 



(81) 



so one obtains 



Ao{k) = 



Ai{k) = 4 



16-.. ]eo^^ 



4^ 
3' 

' Hi 

Hjtk)' 
-"0 



(47r2 + 16C(3)- 39.3139) eg 






(128C(3)- 144.136) ei 



144.5036 -^TT^- 64 C(3) 
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A2{k) = 16eo . 



16 ^Hjtk)^ 



^3(fc) = -el^^ (82) 



with the relation Scq < — W§^ < ^^ which is valid during the inflationary era (which, for us, ends at e(t) — 1/3). 
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